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ABSTRACT 


Statement  of  problem 

This  is  a  general  study  of  the  problem  of  nonequilibrium  drag 
on  an  arbitrary  body  in  steady  adiabatic  subsonic  dissociating  gas 
flow,  viscosity,  heat  conduction  and  diffusion  being  neglected* 


Formulation  of  problem 

Adopting  the  Lagrange's  approach  we  consider  a  moving  material 
element  along  a  streamline.  Besides  the  condition  of  conservation  of 
mass,  the  basic  equations  are 

(1)  thermal  and  caloric  equations  of  state 

(2)  chemical  relaxation  equation 

(3)  general  equation  of  thermodynamic  change 

(4)  first  law  of  thermodynamics 

These  basic  equations  describe  completely  the  thermodynamic  behavior 
of  the  moving  material  element  subjected  to  the  effects  of  dissoci¬ 
ation  and  recombination.  We  linearize  these  equations  and  combine 
them  to  obtain 


a}  dtx  dt*  tL  It  JLt  J  0 


(1) 


which  provides  the  governing  relation  between  the  variations  of  f>' 
and  j'  of  the  moving  material  element  on  a  streamtube.  We  have 
adopted  here  the  following  symbols: 
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f  =  f>  -  f>  s  local  pressure  -  free  stream  pressure, 

l/  =  iT  -  iTe  =  local  specific  volume  -  free  stream 
specific  volume, 

*£  =  relaxation  time  based  on  free  stream  value 

af,  ae  =  frozen  and  equilibrium  speed  of  sound,  respectively 
t  =  time  variable 


General  expression  of  drag  on  a  body 

On  a  streamtube,  we  take  two  points  P  and  Q  denoted  by  z  =  - »o 
and  z  :  too  ,  respectively)  which  are  at  the  positions  upstream  and 
downstream  far  from  the  body  where  all  perturbations  almost  vanish. 

We  show  that  the  body  drag,  D,  can  be  computed  as  follows: 

2  *+«o 

D  * '  -k  // J<r /  fJr  (2) 


2*-o« 

where  jj  denotes  the  integration  with  respect  to  all  stream- 

tubes,  the  streamtube  cross  section  being  c/cT  . 


Results.  Clearly,  to  compute  the  drag  we  must  determine  t>'  > 
along  a  streamtube.  In  the  subsonic  flow  problem  the  function  \f' 
between  P  and  Q  can  be  generally  represented  by  a  Fourier  Series: 

f  OB  )  T  ^  ^ 

IT  =  E  <2*  +  I  cosurn-i  (3) 

where  y*  denotes  the  time  interval  during  which  the  material  element 
moves  from  P  to  Q.  The  real  Fourier  coefficients  an,  bn  can  be  deter¬ 
mined  by  consideration  of  the  dynamic  and  kinematic  conditions  appro¬ 
priate  to  the  problem. 
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(u)„Z)  (  if 


>  0 


(6) 


In  equilibrium  flow,  u)^z  — —  0  /)„**  £>^**0  hence, 

D  =  o.  In  frozen  flow,  — —  *c  .  £>*  =  0  again, 

D  =  o.  In  nonequilibrium  flow,  the  subsonic  drag  must  take  the 
positive  sign.  Thus,  we  have  shown: 

(1)  The  body  will  experience  a  positive  nonequilibrium 
drag  in  subsonic  dissociating  gas  flow. 

(2)  Within  the  accuracy  of  Eq.  (4),  the  nonequilibrium 
drag  is  of  the  second  order  in  perturbation  quantities 
an>  bn* 


iii 


We  show  further  that 


(3)  The  specific  entropy  increases  along  a  streamline 
and  is  of  second  order  in  perturbation  quantities. 

(4)  The  effects  of  entropy  increase  along  the  stream¬ 
line  is  to  contribute  to  the  above  drag  formula  a 
correction  term  which  is  of  the  third  order  in 
perturbation  quantities. 

Remarks 

When  we  consider  the  problem  of  nonequilibrium  supersonic  flow, 
the  drag  would  consist  of  supersonic  wave  drag  and  the  nonequilibrium 
drag.  In  such  a  case,  the  disturbance  propagating  along  the  Mach 
waves  are  damped  by  the  nonequilibrium  effects.  Consequently,  the 
wave  drag  may  decrease  from  the  conventional  value. 
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The  Loss  of  Mechanical  Energy  in  the  Flow 
of  a  Dissociating  Gas 


1 


Introduction 

Clarke^  and  Li  and  Wang^^  recently  discussed  the  problem 

of  nonequilibrium  dissociating  gas  flows  past  slender  bodies.  In 

the  case  of  a  slender  body  with  pointed  ends,  these  authors  obtained 

subsonic  nonequilibrium  drag  formulas  which  are  formally  the  same. 

Clarke^  stated  that  in  the  subsonic  range  a  double-pointed  slender 

body  would  experience  a  positive  nonequilibrium  drag.  However,  by 

applying  the  mean  value  theorem  to  the  integral  expression  in  the 

(2) 

slender  body's  subsonic  nonequilibrium  drag  formula,  Li  and  Wang 
demonstrated  that  the  sign  of  the  integral  may  depend  on  the  slender 
body's  cross  sectional  area  distribution.  Physically  speaking,  be¬ 
cause  of  the  chemical  reaction  in  a  dissociating  gas,  part  of  the 
mechanical  energy  in  the  flow  field  is  transformed  into  internal 
energy  of  the  gas.  Thus,  the  effects  of  chemical  relaxation  are  in 
some  sense  similar  to  the  viscosity  effects  which,  indeed,  also 
transforms  the  flow  mechanical  energy  into  heat.  By  this  physical 
argument,  the  subsonic  drag  on  a  body  placed  in  nonequilibrium  dis¬ 
sociating  flow  should  take  the  positive  sign.  In  the  present  paper, 
we  shall  examine  the  general  problem  of  the  loss  of  mechanical  energy 
in  the  flow  of  a  dissociating  gas.  We  shall  establish  the  physically 
interesting  result  that  a  finite  body  would  experience  a  positive 
drag  in  subsonic  nonequilibrium  dissociating  gas  flows. 


Thermodynamic  Considerations 


We  consider  here  the  stationary  subsonic  flow  of  a  dissociat¬ 
ing  gas  past  an  arbitrary  body.  We  neglect  molecular  transport 
phenomena  such  as  viscosity,  heat  conduction  and  diffusion.  Take 
a  small  material  element  moving  with  the  local  velocity  ^  along  a 
streamtube  of  cross  section  d<T  (Fig.  1).  During  a  small  time 
interval  At  ,  the  moving  material  element  occupies  a  volume 


oliT*  =  ^  *t  d<T 


(1) 


The  work  performed  against  the  environment  by  this  material  element 
is 


<$V=  f  d*r* 


(2) 


where  denotes  the  local  pressure  in  inviscid  flow.  The  first 
law  of  thermodynamics  states  that 


&Q  =  dE  +  W  (3) 

where  denotes  the  heat  received  by  the  material  element  and 

dE  denotes  the  change  in  internal  energy  of  the  material  ele¬ 
ment.  In  adiabatic  flow,  &Q-  0  .  Assuming  local  homogeneity, 


we  define  the  local  density  as  follows: 
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where  *y)  represents  the  mass  in  the  volume  \J~  and  IT  is  specific 
volume  i.e.  volume  per  unit  mass.  The  condition  of  conservation 
of  mass  along  the  streamtube,  viz.,  the  continuity  equation,  gives 

»?  -  f  g  <Lt &t  *•  f0U0<^^7^t  "  rf)c  (5) 

where  the  subscript  o  denotes  quantities  in  the  free  stream. 


From  Eq.  (3)  we  obtain 


dE  +  p  J(miT)  =  0 

Je  +  ftu0at<reAt  f>o(ir=o  (3a> 

#  J 

Let  E*  E/>7  denote  the  internal  energy  per  unit  mass.  Eq.  (3a) 
then  becomes 


o(E  +  pdtT  *  0  (6) 

Consider  the  dissociation  reaction  of  a  symmetrical  diatomic  gas  A2: 

Kf  (7) 

+  Ax  -S'/) /  “E  3 

K*\> 

where  each  A2  molecule  is  assumed  to  be  made  up  from  two  A^  atoms. 
The  species  A3  may  be  either  A2  or  A-^ ,  either  is  assumed  to  be 
equally  effective,  and  Kf  and  are  the  overall  reaction  rate 

constants  for  forward  and  reverse  reactions.  The  dissociating  gas 
is  assumed  to  be  a  mixture  of  perfect  gases.  The  thermal  and 
caloric  equations  of  state  can  be  written  as  follows: 
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P^=  ('+°<)RT  J  /?=  -jfc 


(8) 


E'K  + 


^  Ci/ »)]t  +  <*£) 


(9) 


where  /?  denotes  the  universal  gas  constant,  M2  the  molecular 
weight  of  A2>  D  the  dissociation  energy  per  unit  mas,  o(  the  mass 
fraction  of  the  atomic  species,  and  and  the  specific 

heats  of  the  atoms  and  molecule  respectively.  Simple  calculations 
yield 


ct,  =  [f  +  f  l$>]  R 


C^-  C  ^  ^  “  f  ($)j  R 


(10) 


where  the  translational  and  rotational  energies  are  assumed  the 
fully  excited  classical  values  and  vibrational  energies  are 
assumed  the  equilibrium  values  based  on  simple  harmonic  oscillator 
model  in  quantum  theory,  i.e., 
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kT 


(11) 


h  ,  k  being  Planck's  constant  and  Boltzmann's  constant  respectively 
and  ])  being  vibrational  frequency.  For  simplicity,  Lighthill's 
ideal  dissociating  gas  model^^  can  be  adopted  here,  thus  f(6)*  • 

As  in  Refs.  4  and  5,  the  chemical  relaxation  equation  for  the 
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dissociating  gas  can  be  written  as  follows: 

*  * 

where  Z  and  K  are  defined  as 

#  Mi 

4  f *  0 +*) 


(12) 


(13) 


*  _  Mi.  Kf 

+f  Kv  (1^) 


It  is  worthwhile  to  notice  here  that  Eq.  (12)  is  referred  to  the 
moving  material  element.  The  differentiation  with  respect  to  time 
represented  in  Eq.  (12)  should  be  realized  in  the  Lagrange's  sense. 
If  we  wish  to  have  the  Euler's  derivative,  Eq.  (12)  must  be  inter¬ 
preted  as 


i-*)  -  >'ij 


(12a) 


where  V  denotes  the  local  velocity  vector.  In  the  present  paper, 
will  denote  Lagrange's  derivative. 


General  Expression  of  Drag  on  a  Body 

In  Eq.  (3a),  the  work  performed  by  the  material  element  during 
the  time  interval  dt  is  given  as 


hvt~  p.  ue hB  *t(patir) 


Summing  up  with  respect  to  all  material  elements  along  a  streamtube 
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we  have 


2  *  +  0O 


r 

f.  U,  K  d  J 


2m -o O 


ehr 


Integrating  W]_  expression  with  respect  to  all  the  streamtubes  we 
obtain  the  work  performed  against  the  environment  (in  this  case, 
the  body)  by  the  whole  fluid  during  the  time  interval  A~t  .  Thus 
we  have 

2=  +<* 

lV~  f '.U.djff  pJtrdT.  (15) 

z*  -•« 

The  force  acting  on  the  body  by  the  fluid  along  the  uniform  flow 
direction  is  the  drag,  D.  The  work  performed  against  the  body  dur¬ 
ing  the  time  interval  At  can  be  represented  as 


VJ  =  -DU0At 


(16) 


Comparing  Eqs.  (15)  and  (16)  we  obtain 

Z*-M 

D  =  -P./J7  fJrd'T.  (17) 

In  applying  Eq.  (17)  to  calculate  the  body  drag,  we  must  have  a 
knowledge  of  the  complete  flow  field  around  the  body.  To  ascertain 
the  sign  of  the  drag  integral  we  shall  study  the  sign  of  the  in¬ 
tegral 

Z  m  +  tC 

-f 


J 

Z=  -04 


(18) 
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along  a  streamtube.  In  the  case  of  subsonic  dissociating  gas  flow, 
we  shall  examine  the  equations  of  motion  expressed  in  Lagrange's 
sense.  We  shall  show  that  the  integral  in  Eq.  (18)  always  takes 
the  positive  sign  and  therefore  we  shall  conclude  mathematically 
that  the  nonequilibrium  flow  effects  cause  a  positive  subsonic  drag 
on  the  body. 


Entropy  Change  Along  a  Streamtube 


In  a  nonequilibrium  state,  the  general  equation  of  thermodynamic 
change  in  the  moving  material  element  is  expressible  as^^ 

£>Q  *  Tdt&  (19) 

where  T,  S  denote  temperature  and  entropy,  and  and  denote 
the  chemical  potentials  of  atom  and  molecules,  respectively.  Follow¬ 
ing  Clarke(^),  we  have 


where  b(e  represents  the  equilibrium  atomic  mass  fraction  and  is 
expressible  as  follows: 

<21> 


Combining  Eqs.  (19)  and  (20)  we  have  for  adiabatic  flow 

4  = 


which  gives  the  entropy  change  along  a  streamtube  in  nonequilibrium 


8 


dissociating  gas  flow. 


Linearization  of  Equations 

If  the  body  placed  in  the  flow  is  very  thin,  the  temperature, 
the  pressure,  the  specific  volume  can  be  considered  to  be  slightly 
different  from  the  corresponding  quantities  at  infinity  upstream. 
The  atomic  mass  fraction  o<  is  assumed  to  be  slightly  different 
from  the  local  equilibrium  value  .  Furthermore,  we  assume  cxe 
differs  only  slightly  from  <X0  in  the  uniform  free  stream, 
being  necessarily  an  equilibrium  value.  Then,  we  have 

T  -  T0  +  T 

P  -  F.  +  P' 

u-mir.  +  ir'  (23) 

ci  =  c*e-*-  <X  * 

CX’g  —  o(9  ■+■  Clg 

where  the  subscript  "o"  represents  the  quantities  in  the  uniform 
free  stream  and  the  superscript  "/ "  the  perturbation  quantities. 
From  Eqs.  (6)  and  (9),  we  obtain 


bs«(Cr,-cr,)Jj£*[D  +  (c«-CjT]$Z.  +  p4z~° 


In  obtaining  this  equation,  and  C^-C^  are  assumed  inde¬ 

pendent  of  T,  i.e.,  we  take  f (9)  =  constant.  As  has  been  already 
pointed  out,  this  simplification  corresponds  to  Lighthill's  ideal 
dissociating  gas  model. 
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Substituting  Eq.  (23)  into  Eq.  (24)  and  neglecting  the  squares 
of  perturbation  quantities  we  obtain 

[c^*(c,rc4#+[D+(c,,-c^T.](#+4£)  t  fi#-  o  (25) 


Similarly,  neglecting  the  second  order  terms  of  perturbation  quan¬ 
tities,  we  obtain  from  Eq.  (8) 


fiTZT  +  RT° 


Eq.  (12)  can  be  linearized  in  a  similar  manner  to  yield 

A*.  ~  <^e  +  A*L _ 2£_1  (27) 

dt  7r*7r  ~ 

where  X  represents  the  chemical  relaxation  time  which  is  a  positive 
constant  in  the  present  approximate  calculations^^.  From  Eq.  (22), 
we  obtain 

/ 

ds  „  _ « _ 

d  t~~ZR  «.(,-«?)  ~lt  (28) 

where  only  the  linear  term  in  the  coefficient  of  the  right  hand  side 
term  in  Eq,  (22)  has  been  retained.  Combining  Eqs.  (27)  and  (28)  we 
obtain 


dt  T».  (I-**) 


TT  <*'  +  0  (ex'3) 


This  equation  shows  that  the  entropy  increases  along  the  streamtube 
and  is  of  second  order  in  the  perturbation  quantity  o<  .  By  Eq.  (21), 


we  have 
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=  At-  +  is 

v  JP  K  it  VTS  Jf-TT 

(30) 

By  Eq.  (23),  we  then  have 

ir~  +s-jit 

(31) 

where  \e  and  ^ie  denote  (Tp)s  anc^  (j|s')p  respectively  in  the 

uniform  flow.  Eq.  (31)  shows  that  along  the  streamtube  ^jf*~  > 

and  are  closely  related.  By  Eq.  (29),  we  realize  that 

ij§_  is  of  second  order  in  the  perturbation  quantity  and  can  be 

neglected  in  a  linear  theory.  Therefore,  we  have 

(31a) 

Eliminating  and  from  Eqs.  (25),  (26) 

(31a),  we  obtain  the  following  equation: 

,  (27)  and 

,  J'p  ,  M  c Ip'  J  dV  .  ±L  4tr'  -  * 

L1P  +zdt+"  it'*  *  it  0 

(32) 

“here  L  =  +  C^J  ^ 

(33) 

V  =  LC^+  -  C  ( M- «.)  Rj  />. 

(34) 

M  =  i.  +  1- RT.  {C,,~  ZCrjl 

(35) 

If  the  perturbation  specific  volume  IT  is  given  as  a 

function  of  t 

along  a  streamtube,we  can  find  by  Eq.  (32)  the  perturbation  pressure 
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p  as  a  function  of  ~t  •  We  cannot  have  the  actual  functional  form 
/ 

of  1/  unless  the  equations  of  motion  for  the  flow  problem  are  solved. 


Now,  we  shall  take  two  points  P,  Q  on  a  streamline  (Fig.  2). 
These  points  P,  Q  are  at  the  positions  upstream  and  downstream  far 

from  the  body,  where  all  the  perturbations  almost  vanish.  Between 

/ 

these  two  points,  the  function  iT  in  a  subsonic  flow  problem  can  be 
represented  as  a  Fourier  series 

1  +  I  h„cos2Tn^*  (36) 

n*\  V  n-i  V 

where  denotes  the  time  interval  during  which  the  material  element 
moves  from  P  to  Q.  The  Fourier  coefficients  an  and  bn,  in  principle, 
can  be  determined  for  a  given  flow  problem. 

Solution  of  Eg,  (32)  for  f'=  a„sin  cdnt 

/ 

We  shall  study  the  solution  f>  for  a  Fourier  component 

\r  -  an  sin  cO„t  (37) 


where 


zrn 

V* 


Substituting  Eq,  (37)  into  Eq.  (32)  we  have 

LTFt  'z'ilt  ~ s/n =  o  (38) 
A  solution  of  Eq.  (38)  can  be  written  as 
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(x  s//)  -h  Y cos  od't)  (39) 

where  X,  Y  denote  the  integration  constants  which  can  be  determined 
from  the  following  simultaneous  algebraic  equations 

(r^L)x  +  MY  =  -  t«3„a/ 

(40) 

MX-  (ZtO„L)Y  —  —  Aj 

Solving  Eq.  (40)  with  respect  to  X  and  Y,  we  obtain 

(T^nLf+M3- 

(41) 

y  -  ma i- ln) 

These  results  can  be  discussed  in  the  following  cases:  (1)  equili¬ 
brium  flow  (2)  frozen  flow  and  (3)  nonequilibrium  or  relaxing  flow. 

(1)  Equilibrium  flow.  To Jn — t—  Q 

From  Eq.  (41)  we  have 


X  «  - 


JL 

M 


then  we  obtain  by  Eq.  (39) 


Y  -o 


(42) 


P  =  - 


(43) 


13 


We  easily  find  that 

=  v _ g 

d\r'  ~  <*s'/dt  *  ~  M 


(44) 


Eq.  (22)  shows  that  the  equilibrium  flow  is  isentropic. 
Thus,  Eq.  (44)  can  be  restated  as  follows: 


dhr  dir 


(45) 


where  a£  denotes  the  equilibrium  speed  of  sound  in  the  undis¬ 
turbed  free  stream.  Therefore  we  have 


a 


e 


2. 


(2)  Frozen  flow  Ta)„  1  — ~ 
From  Eq.  (41)  we  have 


(46) 


x— -£ 


Y-  6 


After  the  sane  treatment  as  above,  we  obtain 


-i  N/  1 

n  =  —  i/- 
L 


(47) 


(48) 


where  af  denotes  the  frozen  speed  of  sound  in  the  undis¬ 
turbed  free  stream.  From  Eqs.  (46)  and  (48),  we  have 

l  +  ^-[0**.)RD  +  RT.  (C^- tCrj]  («) 

This  can  be  also  written  as 
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where 


m.  -  x. 

(■w)pf  -  ^[(ViCjT.+dt^Dj 

■o^MKc*+<fc*>"  c*Wi+«.> «] 

h  =  E+-E-  denotes  the  specific  enthalpy  and  suffix  "o" 

denotes  free  stream  quantities.  Eq„  (50)  which  relates 

2  2 

the  quantities  a£  and  ae  has  been  previously  obtained 
by  Clarke(7)  who  also  has  shown  that 


(3)  Nonequilibrium  or  relaxing  flow,  T  =  finite 

Now  we  can  rewrite  the  results  in  Eq.  (41)  as  follows: 


j]  Q/ 

[(  -jf ]  i£* 

T*0*  [t?  "  ']  *1 


X  =  - 


Y-  - 


[(T«>,)*+-^]  IT* 
From  Eq,  (39)  we  obtain 


- J^ZJTZ - J-*,n(^t  +  ‘P)  (53) 
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where 


<P  -t. 


-I  TtJ 


ah 


"  (  k1'1) 


(to)*)1 


(54) 


Eq.  (54)  also  reveals  the  following  facts: 


for 


for 


T<0„  =  finite, 


0  <  cp  <. 


JL 

2. 


T  06„=O  or  T(6a  — •—  oo  cf  =:  o 


Comparing  Eqs.  (37)  and  (53),  we  may  conclude  W  and 
.  / 

P  are  out  of  phase  due  to  relaxation  effects.  In 

/  , 

equilibrium  and  frozen  flows,  however,  \T  and  p  are 
in  phase  and  o.  The  consequence  of  the  phase  shift 
due  to  relaxation  effects,  of  course,  is  to  cause  a 
relaxation  drag  as  can  be  seen  from  application  of  the 
above  results  in  Eq.  (18). 

Subsonic  Drag  due  to  Nonequilibrium  Effects 

In  the  simple  case  discussed  in  the  previous  section,  the 
integral  in  Eq.  (18)  becomes 


z-+- 


(y.f  r')  + 

'Z*-i ©  J 

Z-=-K«  *,“**•  /,  ,  J  .  / 


2»  -oc 


(rj]  jC 


>  0 


(55) 
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In  the  general  case,  vT  should  be  expressed  as  in  Eq.  (36).  The 
corresponding  solution,  p  ,  can  also  be  represented  by  Fourier 
series.  Thus,  the  integral  in  Eq.  (18)  can  be  obtained  as  follows: 
z» +•<> 


'  j 


**1^2. 
I  +  I 


(56) 


where 


•j*  <£  (*>*?)  at  (  ~  0  . 


tOaj/  ^  (qj.r)  4  (v-Q 


0 


(57) 


Thus,  the  integral  in  Eq.  (18)  always  take  the  positive  sign.  By 
Eq.  (17),  the  drag  experienced  by  the  body  in  subsonic  flow  is 


D  -  f.Z-  (/  (At*  £  0  (58) 

In  equilibrium  flow,  oOnX  ~~~~  An  *  sr^  hence,  D  =  0. 

In  frozen  flow,  ^  0  again,  D  =  0.  In 

nonequilibrium  flow,  the  subsonic  drag  must  take  the  positive  sign. 


Discussion 

We  have  studied  the  general  problem  of  subsonic  drag  due  to  a 
nonequilibrium  effect.  In  subsonic  adiabatic  flow,  the  work  performed 
against  the  drag  is  transformed  into  the  internal  energy  of  gas.  Be¬ 
cause  of  the  irreversible  chemical  changes  in  the  gas  there  is  an 
increase  in  the  entropy  of  the  gas.  In  the  present  calculations, 
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the  entropy  change  has  been  neglected  in  the  first  approximation. 

A  higher  approximation  may  be  considered  which  includes  the  effects 
of  entropy  change.  Thus,  in  Eq.  (18),  we  write 

+  Z«+»0 

-J  r.jS-j  fl<r‘ 

Zm  —•o  _ 


Po 


-  -  P.  [  I fi/)  -  irVojf -  f '(f)  r '(/)  +  p'(6)/(o) 


(59) 


/  *  '  V 

where  we  can  take  l/  (  V1  „)  *  *  0 

(Fig.  2)  chosen  to  be  such  that 


irP  - 


by  having  P,  Q 
.  By 


Eq.  (31),  we  have 

7  "  r' 

,  r  L 

0 


(60) 


In  the  linear  theory  in  the  preceding  sections,  the  term  contain¬ 
ing  in  Eq.  (60)  has  been  neglected  because  *  0(«‘  ) 

Within  the  accuracy  of  this  approximation,  we  have  obtained  the 
drag,  the  leading  term  of  which  is  of  the  second  order  in  perturba¬ 
tion  quantities.  If  we  take  account  of  the  increase  in  entropy  of 


18 


the  dissociating  gas  we  may  obtain  a  second  approximate  theory 
which  would  be  correct  to  terms  of  the  third  order  in  perturba¬ 
tion  quantities. 

In  the  conventional  supersonic  flow,  disturbances  propagate 
along  the  Mach  waves.  Mechanical  energy  is  carried  away  to  in¬ 
finity  along  these  waves.  A  part  of  the  conventional  supersonic 
drag,  i.e.,  wave  drag,  must  be  associated  with  the  mechanical 
energy  loss  of  these  wave  disturbances.  When  we  consider  the 
problem  of  nonequilibrium  supersonic  flow,  the  drag  would  consist 
of  the  supersonic  wave  drag  and  the  nonequilibrium  drag.  In  such 
a  case,  the  disturbances  propagating  along  the  Mach  waves  are 
damped  by  the  nonequilibtfium  effects.  As  a  result,  the  wave  drag 
may  decrease  from  the  conventional  value.  Unless  the  flow  problem 
is  completely  solved,  it  would  be  rather  difficult  to  forecase  the 
total  effects  of  nonequilibrium  on  the  supersonic  drag. 
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